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Abstract 

Kv1.3 channel is a delayed rectifier channel abundant in human T lymphocytes. Chronic inflammatory and autoimmune 
disorders lead to the over-expression of Kvl .3 in T cells. To quantitatively study the regulatory mechanism and physiological 
function of Kvl .3 in T cells, it is necessary to have a precise kinetic model of Kvl. 3. In this study, we firstly established a 
kinetic model capable to precisely replicate all the kinetic features for Kvl. 3 channels, and then constructed a T-cell model 
composed of ion channels including Ca^^-release activated calcium (CRAC) channel, intermediate K^ (IK) channel, TASK 
channel and Kvl. 3 channel for quantitatively simulating the changes in membrane potentials and local Ca'^^ signaling 
messengers during activation of T cells. Based on the experimental data from current-clamp recordings, we successfully 
demonstrated that Kvl .3 dominated the membrane potential of T cells to manipulate the Ca^^ influx via CRAC channel. Our 
results revealed that the deficient expression of Kvl .3 channel would cause the less Ca^^ signal, leading to the less efficiency 
in secretion. This was the first successful attempt to simulate membrane potential in non-excitable cells, which laid a solid 
basis for quantitatively studying the regulatory mechanism and physiological role of channels in non-excitable cells. 
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Introduction 

Kvl. 3 and IK (KcaS.l) are two kinds of potassium channels in 
T cells. Kvl. 3 channels are activated upon the depolarization of 
the membrane potential, while IK channels are activated by the 
Calcium ion [1—3]. There are three types of T cells, Naive T cell, 
TcM and Tem cells [3] . These quiescent cells exhibit a similar K"^ 
channel expression pattern with ~300 Kvl. 3 and ~10 IK 
channels per cell. However, Kvl. 3 channels are up regulated to 
~1500 in the activated Tem effectors and IK channels are up 
regulated to ~500 in the activated Tcm effectors [3]. Such high 
expression of Kvl. 3 channels has been reported associated with 
many chronic inflammatory and autoimmune disorders such as 
multiple sclerosis (MS), type 1 diabetes mellitus (TIDM) and 
rheumatoid arthritis (RA) [4—7], and therefore Kvl. 3 channel 
served as a potential therapeutic target for treatment of these 
diseases, which was indicated by the blockers of chemical 
molecules and peptide toxins [8-11]. 

Kvl. 3 channel is a voltage-activated K'*' channel that shows a 
fast activation and slow C-type inactivation and recovery [2,12- 
14]. Comparing with other ion channels, Kvl. 3 is sensitive to 
many pharmacological agents including small organic compounds 
and many peptide toxins such as margatoxin (MgTx), agitoxin-2 
(AgTx2), ShK etc. [15,16]. Particularly, one engineered scorpion 
toxin ADWX- 1 (autoimmune drug from Wenxin group) shows the 
highest affinity to Kv 1 . 3 channels in picomolar range of potency 



[17]. Those properties enable us to identify the Kvl .3 from various 
lymphocyte K"*" currents. 

Besides Kvl. 3, five main types of ion channels have been 
identified at the molecular level in T cells. They are Ca^''"-release 
activated calcium (CRAC) channel, intermediate K''"(IK) channel, 
TASK channel (a two-pore domain potassium (K2p) channel), 
TRPM7 channel and Osmo-activated Cl~ (Cls„(,ii) channel [3] . As 
a major calcium source in T lymphocyte cells, the CRAC channel, 
formed by the STIMl and Orail subunits, leads a Ca^""" influx 
while depleting the endoplasmic reticulum (ER) Ca^"*" store [18- 
20]. Since the intracellular Ca^""" can modulate various important 
physiological functions such as potassium channel gene expression 
and secretion, the Ca influx via CRAC channels may form a 
positive feedback to induce the larger Ca2+ signal [3,18,21]. 

It is well-known that a typical resting effector memory T (T^m) 
cell contains ~300 functional Kvl. 3 channels and ~10 functional 
IK channels on the surface membrane [2,3], which can prevent 
membrane potential from excessively larger depolarization 
[18,19,22-25]. During T cell activation, they are the major 
contributor to maintain the membrane potential that promotes 
Ca ^ influx. A negative membrane potential enhances Ca entry 
by optimizing the electrochemical driving force for Ca^"*" 
movement through CRAC channels [26—28], and inhibit the 
Kvl. 3 and IK potassium channels by their specific inhibitors wiU 
significandy reduce the calcium signaling [8,9]. However, a clear 
pattern in quantitative description on the functional role that 
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Kvl.3 channel plays in regulating the membrane potential and the 
intracellular local and global Ca^""" signaling remains wrapped. 

Kinetic modeling provides a good way for studying and 
predicting the kinetic behavior of ion channels and their function 
in cells. Several kinetic models of Kvl.3 channel have been 
reported previously [12-14]. They did excellent works focused on 
the individual activation, inactivation or recovery characteristics. 
The accuracy of the cell model is based on the comprehensiveness 
of the ion channel models it constitutes. In this study, we establish 
a novel Kvl.3 model capable to precisely describe the whole 
kinetic behavior of K\'1.3, using a software CeL [29]. Based on the 
Hodgkin-Huxley theory, a model cell with appropriate component 
of channels can be used to simulate the firing pattern of action 
potentials in excitable cells [30]. But the H-H model has never 
been used to simulate the membrane potentials in non-excitable 
cells. This is a first attempt to construct a T-cell model, composed 
of several model channels including Kvl .3, CRAC, IK and TASK 
channels, for mimicking the dynamic behavior of membrane 
potentials and intracellular Ca^""" signaling in T cells. Although 
there is no action potential in T lymphocyte cells, it is still 
interesting to know the membrane potential performances after 
stimulation. Combined with the current-clamp experimental data 
with different amount of Kvl.3 channels blocked by ADWX-l 
from T lymphocyte cells, we do rjuantitatively mimic all the 
changes in membrane potential and the corresponding Ca^""" signal 
in the non-excitable model T-cell. Overall, a simulation frame- 
work has been provided for fiarther studying the regulatory 
mechanism of other channels in T lymphocyte cells, ^vhich lays a 
solid basis for both the immunological and ion-channel fields. 

Materials and Methods 

Ethics Statement 

Peripheral venous blood was obtained from healthy volunteers, 
who provided their written informed consent to participate in this 
study. The consent procedure and our research were approved by 
the Ethics Committee of the College of Life Sciences in Wuhan 
University (Permit Number: ECCLS 20120076). 

Cell culture and Transfection 

Full length cDNA for Kvl.3 was subcloned into pcDNA3.1/ 
Zeo (Invitrogen). HEK293 cells were cultured in modified Eagle's 
medium (DMEM, Gibco) supplemented with 10% fetal bovine 
serum (FBS, Gibco) at 37°C incubator with 5% CO2. The day 
before transfection, cells were transfered into a 24-well plate and 
transientiy transfected using lipofectamine 2000 (Invitrogen) 
according to manufacturer's protocol. Recordings were carried 
out in 1-2 days after transfection. 

Preparation of human T lymphocyte cells 

Peripheral venous blood was obtained from healthy volunteers. 
Mononuclear cells were isolated by FicoU-Hypaque density 
gradient. T cells were separated from PBMCs by using FACS 
with PE conjugated CD 3^ antibody and maintained in RPMI 
1640 medium supplemented with 1 mM L-glutamine and 10% 
fetal bovine serum (Gibco) in a humidified, 5% C02 incubator at 
37°C. Recordings were carried out in 24 hours after T lymphocyte 
cells preparation. 

Solutions 

For whole-cell patch recordings, the pipette solution contained 
die following (in mM): 130 KCl, 2.5 MgClj, 10 HEPES, 1 EGTA, 
2 K2ATP (pH 7.4) titrated with KOH, bath solution contained the 
following (in mM): 145 NaCl, 5 KCl, 1 MgClj, 2 CaClj, 10 



HEPES (pH 7.4) titrated with NaOH. All tiie chemicals were 
attained from Sigma. 

Electrophysiology 

Patch pipettes pulled from borosiUcate glass capillaries with 
resistance of 2-3megohms in transfected HEK293 cell experi- 
ments and 6-8megohmsin T lymphocyte cell experiments when 
filled with pipette solution. All experiments were performed using 
an EPC-9 patch-clamp amplifier and corresponding software 
(HEKA, Germany). Currents were typically digitized at 20 kHz 
and filtered at 2.9 kHz (Bessel). The 85% of series resistances were 
electronically compensated. During recording, the corresponding 
solution was puffed onto cells via a puffer pipette containing seven 
solution channels. All experiments were performed at room 
temperature (20-24°C). 

Intracellular Ca^^ Measurements 

T lymphocyte cells were incubated in 100 |J,1 of extracellular 
solution (ECS) containing 140 mM NaCl, 5 mM KCl, 1 mM 
MgCla, 1.8 mM CaCl2, 10 mM Hepes (pH 7.4), 10 mM glucose, 
2 \lM fura-2/AM (Invitrogen), 2 mM Probenecid, 0.05% Pluronic 
F-127 and 0.1% bovine serum albumin at 37°C for 30 min. After 
the incubation, cells were placed in the platform for optical 
imaging. The imaging of intracellular Ca^""" was performed on 
Olympus-IX70 microscope system with a polychromatic light 
source (T.I.L.L. Photonics GmbH, Grafelfing, Germany). The 
excitation wavelengths of fura-2 fluorescence were 340 nm and 
380 nm, coming from the bottom of plate at 1 Hz. AH of the 
experiments were performed at room temperature. 

Data analysis 

Recording data were analyzed with IGOR (Wavemetrics, Lake 
Oswego, OR), Clampfit (Molecular Devices, Inc.) and Sigmaplot 
software (SPSS, Inc.). Unless stated othc-rwise, the data are 
presented as mean ± S.D. The conductance of Kvl.3 channels 
was calculated from G = Ip^^/fV-E/^jg), where Ip^ is the peak 
current, isfo/.j is the reversal potential of Kvl.3. The conductance- 
voltage (G-V) curves of activation were fitted by the Boltzmann 
equation: 

G/Gmax = (l+exp(-(K-F5o)//c))-', (1) 

Where G„ is the maximum conductance, Fjo is the half maximal 
conductance voltage and K is a slope factor. The steady-state 
inactivation curve was fitted to the Boltzmann equation: 

W=/max/(l+exp((K-K5o)/fc)) (2) 

Where I^ax is the maximal peak current, Vso is the half maximal 
availability voltage and the k is a slope factor. 

Mathematical modeling and simulation 

T lymphocyte model cell has been described as the following 
equations [31]: 

Cm * dV/dt = linj + IiK + IcRAC + Ikv13 + ItASK (3) 

where V is the voltage in mV, C„ is the capacity of the cell 
membrane in pF, t is the time in ms, Ii„j is the injection current in 
pA, IfASK is the leak c:urrent in pA, Icrac is the CRAC current in 
pA, I]Q,i3 is the Kvl.3 current in pA and //^is the IK current in 
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pA. In this study, C„=1.5pF and the resting potential 
Vg^, = - 55 mV [32,33]. 
All currents except Ii„j can be written as follows: 

Itask = Gtask * {V —Etask), 

Where the total conductance Gy^iy/; = 0.13 nS [34] and the 
reversal potential Exask- ~ 100 mV [35]; 

IcRAC = GcRAC * {V — EcRAc), 

Where the total conductance Gcrac- 0.09 nS [3] and the reversal 
potential £cft4c = +80 mV [3]; 

JKv\.i = GKyi.i* {V—EKy\.i) 

(4) 

= n * gKvi.i * PoiUV) * {V-EKA.i), 

Where the total conductance Gxvi.s — Variable, n is the number of 
the K^is channel, the a single-channel conductance gK^ij = 1 5 pS 
[36], and the reversal potential Eji:„jg= —13 mV [37], and P„(F) 
comes from the Kyi,3 model. 

hK = GiK * {V-E,K)=n* giK * Po{V,Ca) * (V-Eik), 

where the single-channel conductance giK— 1 1 pS in normal 
Ringer saline [38] , n is the number of the IK channel, the total 
conductance G//, = 0.11 nS, assuming that 10 channels per Tem 
cell [3] and tlu- rc-vcrsal potential -5/^= —75 mV [31], and P„{Ca) 
comes from the P„(Ca)-[Ca^^i curve reported by Grissmer et al. 
[38]. 

Considering that the /A" channel is triggered by a global Ca^""" 
concentration in nM, which usually is tremendously smaller than 
the local Ca^"^ in [39], we thus determined the global Ca 
concentration based on an assumption proportional to the local 
Ca^"*" with a average intracellular Ca^"*" (about 500 nM) after 
CRAC channel opening, which could be measured in T cells [40] . 

To trigger the secretion in T cells, we need to calculate the local 
Ca^^ concentration, depending strictly on net influx and local 
diffusion of Ca^^ near the CRAC channel in T cells. Assuming 
that secretory proteins, i.e. complexes of synaptobrevin, syntaxin 
and SNAP-25 termed SNARE complexes are colocalized to 
individual Ca^"*" sources such as CRAC channel in T cells [41,42], 
we employ a formulation first utilized by Beeler and Reuter [43] to 
convert Ca^"^ influx to a local Ca^"^ concentration: 

d[Ca^+].^ dt = e,rans * IcRAC-ediff * [C(P-*]i (5) 

where [Ca^^i is the intracellular Ca^""" concentration, e,rans is an 
arbitrary transfer coefficient that scales calcium influx to the 
calcium concentration, and eaur is the calcium diffusion coefficient 
[44]. In this study, Ctra™ = 0.006 |xM*pA"'*ms"' and 
eciiiF= 0.003 ms~' [29], which should be able to produce a local 
Ca^"*" concentration high enough (£10 \lM) to trigger secretion 
[45]. Therefore, we can avoid introducing more hypotheses like 
many uncertain diffusion factors, unknown dissociation equilibri- 
um constant and quantity of various buffers and the calcium 
efflux releasing from stores. AU of those factors can be simply 
solved by changing the diffusion and transfer coefficients to reveal 
Ca^"*" properties in T cells. 

The differential equations for the kinetic modeling were solved 
numerically, using a QjMatrix or Five-order Runge-Kutta 



integration method. The fitting procedure is based on a PSO- 
GSS algorithm for direct estimation of rate constants from 
macroscopic currents [29]. The integrating routines were written 
and executed with software CeL (HUST, Wuhan, Hubei, China), 
compiled with the C++ compiler to run under Windows XP [31]. 
KiiK-tic parameters were optimized with CeL as previously 
described [29]. 

Results 

To understand the role Kvl.3 channels play in T cells, it is 
extremely necessary' to establish a kinetic model capable to 
precisely describe its detailed kinetic characteristics under phys- 
iological conditions. After that, we can use this Kvl.3 channel 
model and other four channel models to construct a model T-ceU 
for mimicking the behaviors of membrane potentials and the 
intracellular Ca^"*" signaling and thereby to determine Kvl.3 
channel physiological role in T cells. 

A sequential kinetic model of Kvl.3 channels 

Several non-sequential kinetic models of Kvl.3 had been 
reported for simulating the kinetics of Kvl.3 currents, but they 
focus(^d on the individual activation, inactivation or recovery 
characteristics [2,12-14]. Here we provided a simple model for 
simultaneously simulating the whole kinetics of Kvl.3 currents 
including activation, deactivation, steady-state inactivation and 
recovery (Fig. lA). It is a sequential model of which the activation 
process depends on the forward (activation) rates a and A from the 
closed (Ci) to the open (O) to the inactivated (I) states, the 
deactivation process depends on the backward rates P and B from 
O to C; and the inactivation process depends on the rates r| and (p 
from O to I. AU the rates can be determined by automatically 
fitting the model to the experimental data with the software CeL 
[29]. 

Many studies have confirmed that the heterologously expressed 
Kvl.3 channels possess biophysical and pharmacological signature 
closely resemble that of native Kvl.3 channels in T lymphocytes 
cells [37,46,47]. Their voltage dependence, rates of opening and 
closing, single-channel conductance, inactivation properties, and 
blockade by different compounds including tetraethylammonium 
(TEA) and charybdotoxin (CTX) were indistinguishable from 
native Kvl.3 channels. So we recorded Kvl.3 dynamics from 
transfected HEK293 cells to simplify the modeUng. 

To complete the kinetic model of Kvl.3, experiments were 
performed on Kvl.3-transfected HEK293 cells to acquire the 
activation (Fig. IB), steady-state inactivation (Fig. IC) and 
deacti\ation (Fig. ID), respectively. In Fig. IB, Kvl.3 currents 
(black) exhibited a fast activation and then a slow inactivation, 
evoked by 1 s depolarizing voltage steps ranging from —70 to 
+50 mV in 15 mV increments from a 90 s holding potential of 

— 90 mV to remove the possible inarti\'ation. The forward rates of 
the process C^O^I could be determined by fitting it to the 
activation currents. Here the red lines conferred a perfect fit 
(Fig. IB). The voltage dependence of steady-state inactivation 
currents (black), standing for the availability of channels, was 
obtained by applying a set of conditioning voltages ranging from 

— 80 to — 1 0 mV for 90 s in 5 mV increments and then measured 
at a 150 ms test pulse of 50 mV as indicated (Fig. IC). The 
backward rates of C'^O^I could be determined by fitting it to 
the steady-state inactivation currents. In Fig. IC, all the fits (red) 
were close to the data (black). For deactivation currents (black), 
cells were depolarized to +50 mV for 7 ms from a holding 
potential of —90 mV and then repolarized to various voltages 
between —100 and — 30 mV in 10 mV increments (Fig. ID). 
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Figure 1. Comparison of the kinetic characteristics of Kv1.3 channels between the experimental data and simulations. (A) Kinetic 
model for Kv1.3 channels. The voltage-dependent rates a, fi,r\ and ip are in an exponential form, in which the a, b, c and d are the pre-exponential 
factor in ms"' and the m, n, p and q are the exponential factor in mV; the rates A and B are constants in ms"'; the v is voltage in mV. (B) The 
activation currents of Kvl .3 channels expressed in HEK 293 cells were recorded from a whole-cell patch. Activation was elicited by voltage steps from 
— 70 to 4-50 mV with increment of 15 mV from holding potential at -90 mV applied every 90 s. The voltage protocol is plotted at the bottom. Red 
lines are model simulations (N = 7465). (C) Steady-state inactivation was determined by measuring the peak current at 4-50 mV, following the 60-s 
pre-conditional pulses between —80 mV and —10 mV in increments of 5 mV. The voltage protocol is plotted at the bottom. Red lines are model 
simulations (N=4518). (D) Cells were depolarized to 4-50 mV for 7 ms from a holding -90 mV and then repolarized to various voltages between 
-100 and —30 mV in increments of 10 mV. The voltage protocol is placed at the bottom. Red lines are model simulations (N = 4605). (E) Statistics of 
activation and steady-state inactivation. For activation, V5o= -28.4 ±1.6 mV, s = 4.9±0.5 mV (n = 9); for inactivation, V5o=-43.5± 3.7 mV, 
s = 3.0±0.3 mV (n = 5). (F) Time constants of activation (n = 9) and inactivation (n = 5) were plotted as the function of voltages. (G) Time constants of 
deactivation (n = 10) were plotted as the function of voltages. The red and black lines (or symbols) are for the simulations and experimental data, 
respectively. 

doi:1 0.1 371 /journal.pone.0089975.g001 



Similarly, the backward rate of C^O could be determined by 
fitting it to deactivation currents. AH the fits (red) were firUy 
overlapping to the data (black) (Fig. ID). 

The normalized activation (G- V) curves (black) of Kv 1 . 3 had an 
averaged value of V50 = —28.4 ± 1.6 mV, and the normalized 
steady-state inactivation (availability) curve tested at 50 mV had a 
V50 = — 43.5±3.7 mV (Fig. IE left), consistent with previous work 
on Kvl.3 [37,46,47]. In Fig. IE, both the G-V and steady-state 
curves fi-om simulations (red) were closely coincides with that of 
data (black). Additionally, the time constants of activation (t.J 
(Fig. IF bottom), inactivation (t;) (Fig. IF top), and deactivation 
(Xd) (Fig. IG) were also consistent to the experimental data. 

Based on the perfect fitness to the activation, inactivation, 
deactivation and steady-state inactivation data, all the parameters 
in Kvl.3 model were basically determined. Then we further 
examined whether this model matched the data of recovery from 
inactivation by subtly readjusting the parameters in the Kvl.3 
model. Since the time course of recovery from inactivation is 
mono-exponential, the sequential model with only one inactiva- 
tion state is good enough for the Kvl.3 kinetics. This model 
conferred by and large a rough fit to both the currents recovering 
at —45 and —90 mV. Fig. 2C displayed a set of fractional 



recovery curves arising from a two-pulse protocol (PI = 1 s, 
P2= 150 ms, V=-100, -90, -80 and -45 mV). The time 
constants of recovery were 8.6 s at —45 mV, 15.6 s at —80 mV, 
14.6 s at —90 mV, 14.6 s at —100 mV, approximately consistent 
to the experimental data (Fig. 2D). Interestingly, this result 
indicated that the recovery of Kvl.3 was weakly voltage- 
dependent under physiological condition [14]. So far, a sequential 
kinetic model, well illuminating the whole kinetic behavior of 
Kv 1 . 3 channels, was successfully established. All the parameters in 
model were listed in Table 1. 

Scorpion toxin ADWX-1 is a pore blocker of Kvl.3 
channel without affecting its kinetics 

As we reported previously, the engineered scorpion toxin 
ADWOC-l is a potent specific inhibitor of Kvl.3. It has the highest 
afiinity in ~pM to Kvl.3 channels [17]. In Fig. 3A-B, 30 pM 
ADX'V'X-l blocked the Kvl.3 currents without changing the G-V 
curve (Fig. 3C) and the time constants of activation and 
inactivation (Fig. 3D). This means that the inhibition by 
ADWOC-l only reduces the number of open channels without 
affecting the channel kinetics. In other words, we can explore the 
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Figure 2. Comparison of the recovery of Kv1.3 channels between the experimental data and simulations. (A-B) The representative 
currents of recovery from Inactivatlon were obtained in whole-cell patch mode at -90 mV and —45 mV. Two-pulse protocol: from a holding 
potential of —90 mV, the first pulse 1 (PI) was applied for 1 s at 4-50 mV, and then the second 4-50-mV pulse (P2) was applied for 150 ms after a 
variable interpulse interval at -90 mV in (A) and —45 mV in (B). The detailed recovery currents are shown in the box. The voltage protocol is placed 
at the bottom. Red lines are fits and black ones are data. The channel numbers for simulations are 3788 for —90 mV and 2740 for —45 mV. (C) 
Comparison of the time courses of recovery between the experimental data (black) and the simulations (red) at -100 mV, —90 mV, -80 mV and 
-45 mV. (D) Comparison of time constants of recovery between the experimental data (black) and the simulations (red). The time constants of 
recovery (mean ± SEM) are 8.6±2.7 s (n = 5) at -45 mV, 15.6±1.6 s (n = 3) at -80 mV, 14.6±2.4 s (n = 7) at -90 mV, 14.6±4.7 s (n = 3) at -100 mV, 
respectively. 

doi:1 0.1 371 /journal.pone.0089975.g002 



Kv 1 . 3 role in T cell by simply changing the number of Kv 1 . 3 
channels in applying with different ADWX-1 concentrations. 



Table 1. Parameters of Kv1.3 model. 







Parameter 


Values 


a 


0.448(ms"') 


b 


0.043(ms"') 


c 


0.003(ms"') 


d 


0.00008(ms"') 


A 


280.035(ms"') 


B 


1.648(ms"') 


m 


27.530{mV) 


n 


1 7.528(mV) 


P 


1 74961 (mV) 


q 


1016.330(mV) 



doi:l 0.1 371 /journal.pone.0089975.t001 



Inhibition and characterization of Kv1.3 channels in T 
cells 

Based on the sequential kinetic model of Kvl.3 channels, their 
inhibition features were further investigated in T cells. We first 
examined the voltage-activated currents in T cells in the absence 
and presence of 30 and 300 pM ADWX-1. Our results 
demonstrated that the slow inactivated currents were sensitive to 
ADWX-1 (Fig. 4A). And after applying 300 pM ADWX-1 to T 
cell, there is little outward current left at 50 mV, indicating that 
Kvl.3 current was the major prominent of outward current in this 
case. In current-clamp experiments, we also found that the 
membrane potentials of T cells rapidly went to depolarization 
direction by 15 pA current injection while Kvl.3 channel were 
inhibited by ADWX-1 (Fig. 4B, green and blue), in contrast, the 
membrane potential changed much slowly without ADWX-1 
(Fig. 4B, black). 

To better understand the regulatory mechanism of Kvl.3 
channels in T cells, we constructed a model T-cell, composed of 



PLOS ONE I www.plosone.org 



5 



March 2014 | Volume 9 | Issue 3 | e89975 



Kvl.3 Role in T Cell 



A B 




C D 

1 10000 1 




-50 0 50 -40 -20 0 20 40 60 

Vm (mV) Vm (mV) 

Figure 3. The ADWX-1 toxin blocked the Kv1.3 currents without changing the channel kinetics. (A-B) Kv1.3 activation currents were 
recorded before (black) and after (red) applying with 30 pM ADWX-1. Voltage protocol is placed at the top. (C) The activation G-V curves of Kv1.3 
(n = 6), before (black) and after (red) applying with 30 pM ADWX-1. (D) The time constants of activation and inactivation of Kv1.3 (n = 6), before 
(black) and after (red) applying with 30 pM ADWX-1. 
doi:1 0.1 371 /journal.pone.0089975.g003 



the Kvl.3 channel, Ca^'*"-release activated calcium (CRAC) 
channel, intermediate K""" (IK or KcaS.l) channel and leak 
K''"(TASK) channel [3], to simulate the membrane potentials in a 
T cell, based on the Eq. (3-5) and the open probability Pq of each 
channels (See the Methods and Materials). For more precise, we 
calculated the time- and voltage-dependent total conductance of 
Kvl.3 via the open probability P„ of Kvl.3 kinetic model rather 
than via the formulized term m^ h in classic Hodgkin-Huxley 
model (See Eq. (4)). Essentially, they are the same. According to 
this T-cell model, our simulations revealed that the total 
conductance of Kvl.3 was 4.5 nS (NkvI.3 ~300) for control, 
2.0 nS (Nkvi.3 -133) for 30 nM ADWX-1 and 0.4 nS (Nkvi.3 
-27) for 300 pM ADWX-1 (Fig. 4B). The corresponding 
membrane potential depolarization property in these situations 
depended on the number of Kvl.3 channels, indicating that Kvl.3 
played a dominative role in maintaining the membrane potential 
of T cells. In order to further validate the accuracy of the Kvl.3 
channel model and T cell model, we used one membrane potential 
trace recorded from a T cell as a voltage command to stimulate 
the other T cells, then we obtained a hump-shaped Kvl.3 current, 
which could be well reproduced by the Kvl.3 kinetic model, by 
using the same voltage command (Fig. 4C). This showed that the 
Kvl.3 current started to be activated at about —40 mV and then 
reached to the maximum at ~40 mV, and finally decreased to the 
minimum at more positive voltages due to inactivation. In this 
simulation, we found that the number of Kvl.3 was about 330, 
consistent to the previous report [3] . Note that here we subtracted 



the background current by applying with 1 nM ADWX-1 to get 
the pure Kvl.3 current. 

The membrane potential of T cells is governed by Kvl.3 
channels 

As the T cell is a secretory cell [48] , it is important to know the 
intracellular local and global Ca^"*" signaling. Preferably, it is to 
know the local Ca^""" signaling inside T cells. Since the I-V curve of 
CRAC channel is substantially linear and the membrane potential 
is the only driving force to promote the Ca inflxox [26-28], we 
first calculated the membrane potentials in a model T-ceU 
containing different amount of Kvl.3 channels. Fig. 4D illustrated 
that the changes in membrane potentials of T cells clearly 
depended on the total conductance (or number) of Kvl.3 
channels. Here we chose several typical numbers of Kvl.3 
channel in T cell. For example, the total conductance of 
22.5 nS implied the number of Kvl.3 was 1500 as a single- 
channel conductance of Kvl.3 was 15 pS [36]. An activated T cell 
with 1500 Kvl.3 channels could be considered as an effector 
memory Tem cell [3] . 

Afterwards, we further calculated the corresponding local 
calcium signaling based on Eq. (5). The calculated Ca signaling 
exhibited a Ca^"^ peak at negative potentials, following by a lower 
Ca^^ plateau at positive potentials (Fig. 4E). The maximal local 
Ca^^ concentration around vesicles is up to over 20 |J,M, which 
would be sufficient to induce T-ceU exocytosis [45]. Our results 
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Figure 4. The number of Kv1.3 channels regulates the membrane potential of native T lymphocyte cells. (A) The Kv1.3 currents in a T 
lymphocyte cell, bathed in normal saline, were stimulated by a voltage step from -100 to 50 mV in the whole-cell mode, in the presence of 0 (black), 
30 (blue) and 300 (green) pM ADWX-1 . The inset shows a picture of a T cell. The measured membrane capacitance C^, = 1 .52 pF. The voltage protocol 
is placed at the bottom. (B) A T lymphocyte cell showed the changes in membrane potentials by an injection current of 1 5 pA, in the presence of 0 
(black), 30 (blue) and 300 (green) pM ADWX-1, respectively. A model T cell, composed of Kv1.3, CRAC, IK and TASK channels, conferred the 
simulations (red) corresponding to the changes of membrane potentials resulted from the Kvl .3 conductance of 0.4, 2.0 and 4.5 nS, respectively. The 
resting potential of model cells was —55 mV. (C) The Kvl .3 currents from a native T lymphocytes cell (black) and a model T lymphocytes cell (red) 
were respectively elicited by a voltage command (cyan), conferring a number of Kvl .3 n = 330. Here the background currents were subtracted by 
1 nM ADWX-1. (D) Simulation for membrane potentials in a model T cell stimulated by an injection current of 15 pA, while varying the number of 
Kvl .3 channels. The steady-state values (VsJ of membrane potentials went up to 35.7 mV for 0.4 nS (dark blue), 1 9.8 mV for 2 nS (green), -0.62 mV 
for 4.5 nS (dark yellow), -35.9 mV for 9 nS (pink) and —43.1 mV for 22.5 nS (light blue) at the time of 10 s (gray dotted line), respectively. (E) 
Simulations of intracellular calcium concentrations with the different numbers of Kvl .3 channels. The steady-state calcium concentrations [Ca'^^lsj 
were 3.48 |jlVl for 0.4 nS (dark blue), 4.89 \iM for 2 nS (green), 7.01 i^M for 4.5 nS (dark yellow), 10.47 |iM for 9 nS (pink) and 1 1.48 jiM for 22.5 nS 
(light blue) at the moment of 10 s (gray dotted line), respectively. (F) The time course of the intracellular Ca^* concentrations, stimulated by low 
osmolarity solution labeled in a cyan bar, in a native T lymphocyte cell. The inset shows the Kvl .3 currents recorded from the same T lymphocyte cell. 
doi:1 0.1 371 /journal.pone.0089975.g004 



revealed that the Ca * signaling strongly depended on the 
numbers of Kvl.3 channel. The larger amount of Kvl.3 had 
the higher Ca plateau. In other words, the less Kvl.3, the less 
[Ca^^; for secretion, suggesting that a T cell lacking of Kvl.3 
might lose function. Recently, it was found that the Kvl.3 
knockout mice produced lower incidence and severity of T cell- 
mediated autoimmune encephalomyelitis [49], which was in line 
with our proposed model that the less Kvl.3 channels would 
inhibit the calcium entry into T cells and likely further affected the 
cytokine production. The Kvl.3, IK and CRAC currents in the 
cases of Figure 4 D-4E were shown in Figure 1 . Note that the 
inactivation property of Kvl.3 channel caused a pseudomorph 
paradox of the larger conductance conferring the smaller currents, 
which was actually derived from the smaller open probability due 
to inactivation of Kvl.3 channel. Nevertheless, Kvl.3 currents 
were still much larger than others, indicating that neither CRAC 



nor IK could significantly affect membrane potentials as Kvl -3 
could. This also suggested that the intracellular Ca^*** was 
dominated by Kvl.3 channels. 

The CRAC channel was triggered by IP3-induced depletion of 
the ER Ca^"^ store [18-20]. In Fig. 4F, however, we found that the 
60% lower osmolarity solution induced a global augment of 
intracellular Ca^''" in T cells, consistent to the report by Ross and 
Cahalan [40]. The low osmolarity solution activated a swelling- 
activated CI (Cls„eii) channel by swelling T cell, leading 
depolarization of T cell [3,50,51]. The time course of the global 
Ca concentration showed a Ca peak at beginning, following by 
the higher Ca^"^ plateau, comparable with that as predicted by the 
model T cell (Fig. 4E). Obviously, the slow effect of the regulatory 
volume decrease (RVD) could not be the major reason for the 
Ca^"^ decay in seconds, because it occurred in minutes. Addition- 
ally, it was not clear whether the low osmolarity can directly 
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Figure 5. A putative activation mechanism of T cells. Motile polarized T cells naturally use their T-cell receptor (TCR) to come into contact with 
antigen(Ag)-presenting cell (ARC) via peptide-major histocompatibility complex (MHC); T cells tagged activate the smaller amount of Ca^* influx via 
CRAC channel, probably evoked by IP3-induced depletion of Ca^* stores, to slightly increase the expression of Kv1.3; more Ca'^* influx via CRAC 
channel enters into T cell due to membrane repolarization by augment of Kv1.3; The higher Ca'^* in T cell triggers the IL-2 promoter production 
program. 

doi:1 0.1 371 /journal.pone.0089975.g005 



trigger the opening of tlie CRAC channel or stretch-activated 
cation channel with a reversal potential V„v~0 mV. For 
simulation, there was actually no essential difference between the 
open of CRAC and stretch-activated cation channels at about 
—40 mV (<0 mV) in applying with the lower osmolarity solution 
[3], because the calcium influx derived from either the CRAC or 
cation channel was governed by Kvl.3 channels. 

Discussion 

The effector memory Tem cell, compos ed ofKvl.3'''*=''Kca3.l'°" 
has been confirmed for several human autoimmune diseases, such 
as multiple sclerosis (MS), type 1 diabetes mellitus (TIDM) and 
rheumatoid arthritis (RA) [5-7]. The Kvl.3 voltage-gated potassi- 
um channel regulates membrane potential and calcium signaling in 
Tem cells that are key mediators of autoimmune diseases, a 
thoroughly studying the physiological role of Kv 1 . 3 channel plays in 
T cell not only helps to understand the working mechanism of T 
cell, but also provide a crucial clue for treatment of diseases. 

In this study, we systematically studied the basic characteristics 
of Kvl.3 channel and constructed a sequential kinetic model 
capable to precisely replicate all the kinetics of this channel. 
Moreover, we examined the Kvl.3 currents in T lymphocyte cell 
by a specific inhibitor ADWX-1 toxin, consistent to that of 
expressed in HEK293. To explore the physiological role of Kvl.3 
channel, we further built a T model cell, which was composed of 
four basic channels prominently in T lymphocyte cell. Finally, we 
demonstrated that the Kvl.3 channels could dominate the 
membrane potentials of T lymphocyte cell, and that the CRAC 
channels manipulated the intracellular Ca^^ signaling along with 
changes of membrane potentials. 

Using the model T-ceU composed of four channels, we have 
examined the mechanism of Kvl.3 channel regulating T 
lymphocyte cells. Actually, a couple of channels, i.e., TRPM7, 
Osmo-activated CI (Cl,„,,n), and TRPC3/6 were excluded from 
this T-cell model. The TRPM7 (or a Mg2"*"-inhibitor-Ca2+ (MIC) 
permeable channel) is a stretch- and swelUng-activated cation 
channel, which can be activated by PIP2 and inhibited by Mg^'*" 
[52,53]. The Cls„eU channel conducts an outwardly rectifying 
chloride current (or volume-regulated anion current), leading an 
anion and osmolyte efflux that is activated by cell swelling [54] . In 
our experiments, TRPM7 current was blocked by intracellular 
2.5 mM Mg^"^. Besides CRAC channel, TRPC3/6 channels 
provide another Ca^"*" entry pathway, but it seems that these two 



channels are littie expressed in T cells, the amplitude of TRPC3 
inward current is only around 2 pA/pF (or 1-3 p A/ cell) at 
— 50 mV [55], which could be neglected in our model, so the 
model T-ceU composed of four basic channels is thus good enough 
for this work. 

During the whole depolarization process of membrane poten- 
tials, CRAC contributes much less current than that of Kvl.3, due 
to the small total conductance of 0.09 nS compared with ~4.5 nS 
of Kvl.3. And the membrane potential was much less sensitive to 
the CRAC channel conductance in simulations, suggesting that 
Kvl.3 played the key role in dominating membrane potentials in 
T cells. While CRAC channel was the major Ca^""" source in T 
cells, which provides a continuous calcium influx to elevate the 
intracellular Ca^""" level for T-cell exocytosis. With TASK channels 
together, Kvl.3 channels also play a role in maintaining the rest 
potential. 

In this study, we identified the physiological role of Kvl.3 
channels, which further provided us an opportunity to better 
understand the activation mechanism of T cells. How can a 
CRAC channel in a T cell be activated following antigen 
stimulation? In Fig. 5, a cartoon shows that an activated T cell 
increases the surface expression of Kvl.3, to negative-shift the 
membrane potential to increase the intracellular Ca^'*' via CRAC 
channels, and finally to start the IL-2 (interleukin 2) promoter 
production program [3]. This process is more like a positive 
feedback [3,18], possibly terminated by endocytosis due to 
internationalization of Kvl.3 channels. This will decrease the 
number of Kvl.3 on plasma membrane to lead T cell gradually 
return to its resting state again. 

In this study, we provided a precise model of Kvl.3 and a novel 
non-excitable model cell for quantitatively investigating the 
regulatory mechanism of T cell, which laid a solid basis for 
further studying other channel such as IK, TASK roles in T cell 
and the events of membrane potentials, calcium signaling and 
secretion in other non-excitable cells such as Tqm; B and sperm 
cells. 
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